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Generalized vector form factors of the pion, related to the moments of the generalized parton dis- 
tribution functions, are evaluated in the Nambu-Jona-Lasinio model with the Pauli-Villars regular- 
ization. The lowest moments (the electromagnetic and the gravitational form factors) are compared 
to recent lattice data, with fair agreement. Predictions for higher-order moments are also made. 
Relevant features of the generalized form factors in the chiral quark models are highlighted and the 
role of the QCD evolution for the higher-order GFFs is stressed. 
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I. INTRODUCTION 

Generalized Parton Distributions (GPDs) are very in- 
teresting physical objects, as they encode in a natural 
way the rich information on the internal structure of 
hadrons (for extensive reviews see e.g. 
and references therein) . Moments of the GPDs in the X 
variable form polynomials in the ^ variable (see the fol- 
lowing), with the t-dependent coefficients known as the 
generalized form factors (GFFs). Experimentally, the 
GPDs, especially for the pion, are rather elusive quan- 
tities, as they show up in difficult to measure hard exclu- 
sive processes such as Deeply Virtual Compton Scatter- 
ing (DVCS) or hard electroproduction of mesons (HMP). 
On the other hand, the GFFs are accessible to lattice 
QCD measurements. In particular, the lowest-order vec- 
tor and tensor [l3| GFFs of the pion have recently 
been determined from fuU-QCD lattice calculations and 
more similar results are announced to appear in the fu- 
ture. 

In this paper we evaluate the vector GFFs of the pion 
in a chiral quark model, namely the Nambu-Jona-Lasinio 
(NJL) model with the Pauli-Villars (PV) regularization, 
using the techniques described in Ref. |11| . Importantly, 
the results for the GPDs obtained in this model obey all 
formal requirements, such as normalization, proper sup- 
port, polynomiality and positivity constraints. At the 
same time, the obtained expressions have a non-trivial 
form which is not of a factorizable in the t-variable. We 
compare the lowest form factors, namely the electromag- 
netic and the gravitational form factor, to the recent lat- 
tice data ^ and find good qualitative agreement. Pre- 
dictions for the higher-order form factors, which will be 
be confronted to data when more lattice results appear. 
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are also shown. A careful analysis of the role of the QCD 
evolution is presented, which is essential for higher-order 
GFFs. 



II. GENERALIZED PARTON DISTRIBUTIONS 
AND GENERALIZED FORM FACTORS 

We begin with the description of kinematics of GPDs. 
The assignment of the momenta in the considered process 
is depicted in Fig. [T] The standard notation is 

= ml, = -2p- 9 = t, 

— 0, p ■ n — 1, q ■ n — —<^. (1) 

The null vector n defines the light cone. Throughout this 
paper we work for simplicity in the strict chiral limit with 

= 0, (2) 

although an extension to the physical pion mass is 
straightforward in the NJL model. The two isospin pro- 
jections of the GPDs of the pion are defined through 
the matrix elements of bilinears of quark fields displaced 
along the light cone, namely 

Sabn'-'>ix,Ct) - J ^e"^'^^" X (3) 

tes,,H'=\x,C,t) = J ^e"f^^" X (4) 

where z is on the light cone, and a and b are the 
pion isospin indices. At the scale pertaining to the chi- 
ral quark models the gluons are integrated out, hence 

In chiral quark models at the leading- iVc level the cal- 
culation of the GPDs proceeds according to the one- 
loop diagrams of Fig. [TJ Extensive details of the quark- 
model evaluation are given in [ll| . The pion GPDs were 
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FIG. 1: The direct (a), crossed (b), and contact (c) Feynman 
diagrams for the quark-model evaluation of the GPD of the 
pion. The contact contribution is responsible for the D-term. 



also analyzed in other variants of chiral quark models in 
Refs. [if M, Q [H, [H, In this paper we use the 
non-linear NJL model with the PV regularization in the 
twice-subtracted version proposed in Ref. p^ . where an 
observable O is regularized according to the prescription 



acg = O(0)-O(A2) + A 



dA2 ■ 



(5) 



8^ 



(7) 



^2^/4 (M2-KA2)-tlog 



-y4(M2+A2)_t_/3f 
^4(M2+A2)-t+v^ 



V 



The condition lim(^_oc Fnjl(0 = is satisfied due to 
Eq. ®. 

In the so-called symmetric notation for the GPDs, 
more convenient for our study, one introduces 



c 



X = 



x-C/2 



2-C' " l-C/2' 
where < ^ < 1 and —1 < X < 1. Then one defines 



(8) 



(9) 



which exhibit the reflection properties about the X = 
point, 

if^=0(X,C,0 = -H'=''{-X,^,t), 
H'=\X,^,t) = H'-\^X,^,t). (10) 

In addition, for X > one has 

F^=°-l(X,0,0)=q(X), 

relating the distributions to the the pion's forward diag- 
onal parton distribution function (PDF), q{X). 

The polynomiality conditions [l], 0] state that the mo- 
ments of the GPDs can be written as 

.1 3 

/ dX X'^ H'=\X, Lt)^2j2 A2,+i4m'\ 
dX X^^+^ H'=°{X,^,t) = 2j2A2j+2At)e\ (11) 

•^-1 1=0 



The quantity A is the PV regulator. In what follows we 
take M = 280 MeV for the constituent quark mass and 
A = 871 MeV, which yields / = 93.3 MeV for the pion 
decay constant [l^ according the the formula 

Calculations in other variants of the chiral quark models, 
in particular in the Spectral Quark Model [l^ [2^ , will 
be presented elsewhere. They are qualitatively similar to 
the results presented here. 

The pion electromagnetic form factor in the NJL model 



where A2j+i,i{t) are the generalized form factors, de- 
pending on = 0, 1, . . . and i. The polynomiality prop- 
erty follows from very basic field-theoretic assumptions 
such as the Lorentz invariance, time reversal, and her- 
miticity, hence is automatically satisfied in approaches 
that obey these requirements. In our approach polynomi- 
ality is manifest from the use of the double distributions 




The notation and normalization factors in Eq. pip are 
adjusted in order to agree with the conventions of Ref. 0], 
except for the subscript i which in our case labels the 
powers of and not ^. We note that for the isovector 
(non-singlet) GPD only even, and for the isoscalar (sin- 
glet) GPD only odd moments are non-zero. For the few 
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FIG. 2: The vector and gravitational form factors of the pion 
in the NJL model with PV regularization in the chiral limit 
(lines) compared to the lattice data from Ref. ^ (statistical 
errors of the data not displayed). 



lowest values of j we have explicitly 



fdX X H'-''{X,^,t) = A2„{t) + A2i{t)e 
Jo 



(12) 



(13) 



fdXX^H'=^ {X, t) = Aso (t) + A31 {t)e 
Jo 



dX X^H'=\X, e, t) = A4o(t) + A4i(t)e' + Ai2{t)e- 



In Eq. we have introduced is the electromagnetic 
form factor Fv{t) = Ai^{t), while 9i{t) = 2A20 and 
6*2 (i) = —2^21 in Eq. are the gravitational form fac- 
tors of the pion, discussed in more detail in In the 
chiral limit these form factors satisfy the low energy the- 
orem 01(0) = ^2(0) '^I'l. In our quark model calculation 
in the chiral limit 



eiit) = 92it) = e{t), 



(14) 



hence, consequently, A2o{t) = —A2i{t). The sum rule 
(fT^ expresses the electric charge conservation, while 
is responsible for the momentum sum rule in the deep 
inelastic scattering. 

Formulas (|lip are equivalent to the definition 

(7r+(p')| w(0)7^''iS'^MS'"" .. .il?^^^u(0) \tt+{p)) = 
2p{MpPi pw}Aj+i.o(t) + (15) 
j 

2 J2 Q^^'q''' ■ ■ ■ q'-'P'-'^' ■ ■ ■ P^-^ A, + i,(,+i)/2(t), 
1=1 

odd 

where P = {p + p')/2, D = \{D-D), and {. . .} denotes 
symmetrization and subtraction of traces for each pair of 



indices. Equivalence of Eq. (I15p and ((TT|) is easily seen 
by contracting fTS]) with the null vectors n^^ . . . n'^J and 
applying the definitions ([T]). 



III. RESULTS AT THE QUARK-MODEL SCALE 

The calculation of the GPDs made in the NJL model 
according to the diagrams of Fig. [1] is straightforward [ll| 
and is most efficiently done via the double distributions. 
Then the the GFFs are extracted from the full GPDs by 
evaluating the moments (|lip . The obtained expressions 
are rather lengthy, hence we do not present them here. 
They have the form similar to Eq. ([7]), involving logs and 
rational functions in the t variable. 

Our determination of GPDs and GFFs corresponds to 
the quark model scale Qo, where matching to QCD is 
made. The reader is referred to Ref. [Tlj , where the issue 
is discussed in detail. The value of Qo may be estimated 
be performing the evolution of the parton distribution 
functions (PDF) or the parton distribution amplitude 
(PDA) to higher scales and comparing the results to the 
available data. It turns out to be low, Qo — 320 MeV. 

The results for the electromagnetic and gravitational 
form factors are independent of the scale. They are 
shown in Fig. O The dots are the data points from lat- 
tice calculations of Ref. [9|. One has to bare in mind 
that the fuU-QCD lattice calculations need to be extrap- 
olated to the chiral limit. Also, our model incorporates 
only the leading- A^c contributions. Nevertheless, we note 
a rather remarkable qualitative agreement, in particular 
the feature of a much slower decay of the gravitational 
form factor compared to the electromagnetic form factor 
occurs both in the model and the data. Of course, for the 
electromagnetic form factor accurate experimental data 
could could be used for the comparison. However, the 
focus of this work is on higher-order GFFs, where the 
information comes from the lattices. 

In Fig. [3] we show the higher level GFFs, A^^i and Ai^i 
obtained at the quark model scale. Since in chiral quark 
models in the chiral limit one has at t = [ll|, [H, l2^ 



H'=\X,^,0) ^9{l-X^) 
H'=°{X,^,Q) = 

0((1 - A)(A - 0) - e((A + 1)(-^ - A)), 



(16) 



it follows from the definition (fTTj) that at the quark-model 
scale 



^2j-Hl,z(0) 



(0) 



2J+1 fori = 
otherwise 

for i = 

2J+2 fori=j + l 
otherwise 



2 J +2 



(quark — model scale) (17) 



This behavior is seen in Fig. [31 We note that the model 
GFFs go to zero very slowly at large —t. Another prop- 
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erty follows from the fact that in the considered model 
H'=°{X, l,t)=0 for any value of t. Then Eq. lU) yields 



i=0 



(18) 



This feature can be seen in the lower panel of Fig. [3l 



IV. QCD EVOLUTION 

As already mentioned, the crucial role of the QCD evo- 
lution in chiral quark model calculations has been dis- 
cussed in Ref. [11] . We carry the leading-order DGLAP- 
ERBL evolution from the quark-model scale 



Qo = 313 MeV 



(19) 



to the scale of the lattice calculation of Ref. |9|. This 
scale can be inferred from the value 2A2o{t) = 0.63 in 
, which is reproduced in our calculation when we evolve 
the isoscalar GPD to the scale = o.71 GeV^. Thus 
we estimate the lattice scale as 



Q = 843 MeV. 



(20) 



We use the method and code described in [23| to evolve 
the GPDs taken at several selected values of ^. From this 
one may disentangle the coefficients of the powers of ^, 
namely the GFFs, at the lattice scale Q. The results of 
this procedure are presented in Fig. [H We note a sizable 
change, both in the value at t = and in shape, compared 
to the behavior of Fig. [31 which shows that in general the 
GFFs do evolve with the scale, except for the protected 
form factors as those in Fig. [21 which are invariants of 
the evolution. 

For the evolution scale ^ oo, the GPDs tend to 
their asymptotic forms located entirely in the ERBL re- 
gion \X\ < ^. Explicitly, we have in this limit pl| 



H 



1=1 



|(i-flM<) 



(21) 



7=0 



H 



XH 



15 



N 



X ( X^" 



f 



X^ 



where Cp = {N^ - l)/{2Nc) and TV/ = 3 is the num- 
ber of active flavors. The proportionality factors reflect 
the normalization at the initial quark- model scale Qq, as 
the charge- and momentum-conservation sum rules are 
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invariants of the evolution, 



J ^dXH'-\X,tt,Q^)^2Fv{t), (22) 
J ^ dX {XH'-\X, t, Q") + XHg{X, e, t, Q')) 



in accordance to Eq. ()12|13|1 . 

Evaluation of moments in Eq. pT|) yields for ^ oo 



[ otherwise 

{4CF + Af/ 2[4i(i+4) + 15]^(^) for Z = j 
-^2i+2 J (i) for I = j + 1 

otherwise 

(Q' oo) (23) 
with the gluon form factors defined as 
.1 i+i 

/ dXX^^+^Hg{X,U,Q^) = 2j2A%+2At)^''-i2^) 

We note a striking difference of the asymptotic form fac- 
tors compared to the form factors at the quark model 
scale shown in Fig. [3] For the isovector case only the 
highest form factor, for i = j, is non-zero, and for the 
isoscalar case only the two highest moments, for i = j + 1 
and i ^ j are non-zero. In contract, at the quark- 
model scale all GFFs are present. Thus only and 
A20 = —A21 are invariants of the QCD evolution. The 
asymptotic gluon form factors in Eq. (|23p are related 
to the isoscalar quark form factors in a simple manner. 
Asymptotically, all GFFs become proportional to Fv{t) 
or 9{t) in the isovector and isoscalar channels, respec- 
tively. The universality for the gluon form factors follows 
from the fact that we start from the quark-model initial 
condition, which carries no gluons, which are then built 
in the process of evolution. Asymptotically, for A^'c = 3 
and Nf = 3 the quark to gluon ratio in the momentum 
sum rule equals 9/16. 



V. SUMMARY 

In summary, we have computed the generalized vec- 
tor form factors of the pion in the Nambu-Jona-Lasinio 
model with the Pauli-Villars regularization. We have 



proceeded through the X-moments of the generalized 
parton distribution functions, evolved from the quark- 
model scale to the scale of the lattice calculations. The 
model GPDs exhibit no factorization in the t-variable. 
Comparison to the lattice results for the electromagnetic 
and gravitational form factors of the the pion have been 
made, with proper agreement. The QCD evolution has 
been carried out and its role for the higher-order gener- 
alized form factors has been discussed. Our predictions 
for higher-order GFFs may be compared to lattice results 
when these become available. 



VI. MANOJ 

Since this volume is devoted to the memory of Manoj 
K. Banerjee, our unforgettable teacher and friend, let me 
finish with an anecdote. When I started my graduate 
work, Manoj gave me a code written together with Mike 
Birse who was then a postdoc at Maryland. The code 
was used to solve the chiral soliton model with valence 
quarks (the Birse-Banerjee model [l^, H^), which was 
a major achievement. My assi gnm ent was to introduce 
vector mesons to the model j26l.l27l| . When reading down 
the FORTRAN lines I noticed that the encoded radial 
differential equation for the pion had a seriously-looking 
mistake: in one of the terms instead of p/r**2 there 
was p/r*2. Omission of one asterisk changed the square 
into the multiplication by 2, a potentially devastating 
error. When I showed this to my advisor, he rushed 
to the computer and ran the corrected code. The first 
output quantity was the soliton mass, which changed by 
a tiny amount at the relative level of 10"** or so. Other 
observables were also practically unaffected ... - "You 
see, young man, good physics is immune to such silly 
mistakes as omission of an asterisk!" 

Let Manoj 's kindness, brilliance, enthusiasm, and con- 
fidence be with us! 
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